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The signature of local cosmic structures on the ultra-high energy cosmic 
ray anisotropics 

A. Cuoco''* 

'^Dipartimento di Scienze Fisiche, Universita di Napoli "Federico 11" , and INFN-Sezione di Napoli, 
Complesso Universitario di Monte Sant'Angelo, Via Cintia, 1-80126 Napoli, Italy 

Current experiments collecting high statistics in ultra-high energy cosmic rays (UHECRs) are opening a new 
window on the universe facing the possibility to perform UHECR astronomy. Here we discuss a large scale 
structure (LSS) model for the UHECR origin for which we evaluate the expected large scale anisotropy in the 
UHECR arrival distribution. Employing the IRAS PSCz catalogue as tracer of the LSS, we derive the minimum 
statistics needed to reject or assess the correlation of the UHECRs with the baryonic distribution in the universe, 
in particular providing a forecast for the Auger experiment. 



1. Introduction 

At energies above a few x 10^^ eV, which we will 
refer to as the ultra-high energy (UHE) regime, 
protons propagating in the Galaxy retain most 
of their initial direction. Provided that the ex- 
tragalactic magnetic field (EGMF) is negligible, 
UHE protons will therefore allow us to probe into 
the nature and properties of their cosmic sources. 
However, due to quite steep CR power spec- 
trum, UHECRs are extremely rare (a few parti- 
cles km~^ century"^) and their detection calls for 
the prolonged use of instruments with huge col- 
lecting areas. One further constraint arises from 
an effect first pointed out by Greisen, Zatsepin 
and Kuzmin [1,2] and since then known as GZK 
effect: at energies E ^ 5 x 10^^ eV the opac- 
ity of the interstellar space to protons drastically 
increases due to the photo-meson interaction pro- 
cess p-l-7cMB — > 7r°(+) +p{n) which takes place on 
cosmic microwave background (CMB) photons. 
In other words, unless the sources are located 
within a sphere with radius of O(IOO) Mpc, the 
proton flux at E }^ 5 x 10^^ eV should be greatly 
suppressed. However, due to the very limited 
statistics available in the UHE regime, the ex- 
perimental detection of the GZK effect has not 
yet been firmly established. 

A detailed knowledge of the UHECR frame- 
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work is still missing. Only recently, for instance, 
magnetic fields were included in simulations of 
large scale structures (LSS) [3,4]. Qualitatively 
the simulations agree in finding that EGMFs are 
mainly localized in galaxy clusters and filaments, 
while voids should contain only primordial fields. 
However, the conclusions of Refs. [3] and [4] are 
quantitatively rather different and it is at present 
unclear whether deflections in extragalactic mag- 
netic fields will prevent astronomy even with UHE 
protons or not. Another large source of uncer- 
tainty is our ignorance on the chemical compo- 
sition of UHECRs, mainly due to the need to 
extrapolate for decades in energy the models of 
hadronic interactions, though important progress 
are expected from new high quality data and de- 
convolution techniques (see e.g. [5]). Future ac- 
celerator measurements of hadronic cross sections 
in higher energy ranges will also ameliorate the 
situation, but this will take several years at least. 

From now on, therefore, we shall work under 
the assumptions that UHE astronomy is possible, 
namely: i) proton primaries; ii) EGMF negligibly 
small; iii) extragalactic astrophysical sources are 
responsible for UHECRs acceleration. A possi- 
bility favoring these hypotheses is that relatively 
few, powerful nearby sources are responsible for 
the UHECRs, and the small scale clustering ob- 
served by AGASA [6] may be a hint in this direc- 
tion. However, the above quoted clustering has 
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not yet been confirmed by other experiments with 
comparable or larger statistics [7,8], and prob- 
ably a final answer will come when the Pierre 
Auger Observatory [9] will have collected enough 
data. Independently on the observation of small- 
scale clustering, one could still look for large scale 
anisotropies in the data, eventually correlating 
with some known configuration of astrophysical 
source candidates. In this context, the most nat- 
ural scenario to be tested is that UHECRs cor- 
relate with the luminous matter in the "local" 
universe. This is particularly expected for candi- 
dates like gamma ray bursts (hosted more likely 
in star formation regions) or colliding galaxies, 
but it's also a sufficiently generic hypothesis to 
deserve an interest of its own. 

To this aim we use the IRAS PSCz [10] astro- 
nomical catalogue as tracer of Large Scale Struc- 
tures from which, once that the related selection 
effects have been taken into account, we derive 
the related pattern of anisotropics; wc then as- 
sess the minimum statistics needed to detect the 
model against the isotropic null hypothesis, in 
particular providing a forecast for the Auger ex- 
periment. Previous attempts to address a similar 
issue can be found in [11-13]. Further details of 
the analysis summarized here can be found in [14]. 

2. Data and Formalism 

2.1. The Catalogue 

The IRAS PSCz catalogue [10] contains about 
15.000 galaxies and related redshifts with a well 
understood completeness function down to ^ ^ 
0.1 (i.e. down to a redshift which is comparable 
to the attenuation length introduced by the GZK 
effect) and a sky coverage of about 84%. The 
incomplete sky coverage is mainly due to the so 
called zone of avoidance centered on the galactic 
plane and caused by galactic extinction and to a 
few, narrow stripes which were not observed with 
enough sensitivity by the IRAS satellite (see [10]), 
These regions are excluded from our analysis with 
the use of the binary mask available with the 
PSCz catalogue itself. The mask /i is defined so 
that /Lt(f^) = 0, 1 if the direction H, falls respec- 
tively inside or outside the blind region. 

To correctly employ the catalogue, selection 



and discreteness systematic effects have to be con- 
sidered. For the flux selection correction the rele- 
vant quantity to be taken into account is the frac- 
tion of galaxies actually observed at the various 
redshifts, a quantity also known as the redshift 
selection function (f>{z); given 4>{z), the quantity 
n{z)/<j){z) represents the experimental distribu- 
tion corrected for the selection effects, which must 
be used in the computations. The discreteness ef- 
fects (or "shot noise") are, instead, minimized us- 
ing the catalogue only until a maximum redshift 
of z = 0.06 (corresponding to 180 h~^Mpc), a fair 
compromise for which we have still good statistics 
while keeping the intrinsic statistical fluctuations 
under control. In any case, due to GZK effect in 
the energy range E >5xl0^^ eV, the contribution 
from sources beyond z ~ 0.06 is sub-dominant, 
thus allowing to assume for the objects beyond 
z = 0.06 an effective isotropic source contribu- 
tion. A detailed discussion of the various cata- 
logue issues can be found in [14] and references 
therein. 

2.2. Analysis 

The goal of our analysis is to obtain the under- 
lying probability distribution fhSsi^- E) to have 
a UHECR with energy higher than E from the 
direction Cl. 

For simplicity here we shall assume that each 
source of our catalogue has the same probabil- 
ity to emit a UHECR, according to a power law 
spectrum at the source E~^; however, both as- 
sumpitions can be easily generalized within the 
same formalism [13]. 

A suitable expression for the distribution 
/Lss(f^)-E) can be easily obtained as 

/Lss(fi,i;eut) oc Y: J , T.'^.^i r ^"^di; 

^ 4>{zk) 47rd|^(zfe) Je4e,^,.z,,) 
= Y.fLSs{k)6{h-nk), (1) 

k 

that can be effectively seen as if at every 
source k of the catalogue it is assigned a weight 
/Lss(fc) that takes into account geometrical ef- 
fects through the luminosity distance {dj^'^), se- 
lection effects {4>~^), and physics of energy losses 
through the integral in dE; in this "GZK inte- 
gral" the upper limit of integration is taken to be 
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infinite, though the resuh is practicaUy indepen- 
dent on the upper cut used provided it is much 
larger than 10^° eV; the energy losses are taken 
into account through the code described in [14] 
and parameterized in the function Ei{Ef, z), the 
"initial energy function" giving the energy Ei of 
particle that should be injected at a redshift z to 
reach the Earth with an energy Ej. 

We choose Ecnt = 5 x 10^^ eV that results in 
a fair compromise between the intensity of the 
anisotropics (that increases with energy) and the 
achievable statistics; for this £^cut the isotropic 
contribution to the flux is sub-dominant; how- 
ever we can take it exactly into account and the 
weight of the isotropic part is given by Wiso cx: 

Finally, to graphically represent the result, the 
spike-like map of Eq. (1) is effectively smoothed 
through a gaussian filter with beamwidth u = 3° . 

Given the extremely poor statistics available 
in UHECR astrophysics, we limit ourselves to 
address the basic issue of determining the mini- 
mum number of events needed to significatively 
reject "the null hypothesis". To this purpose, 
it is well known that a x^-test is an extremely 
good estimator and has no ambiguity due to the 
2-dimensional nature of the problem respect to 
the K-S test or the Smirnov-Cramer-von Mises 
test. A criterion guiding in the choice of the 
bin size is the following: with N UHECRs events 
available and M bins, one would expect 0{N/M) 
events per bin; to allow a reliable application of 
the x^-test, one has to impose N/M > 10. Each 
cell should then cover at least a solid angle of 
Am 10 X Atot/-^^, Atot being the soHd angle 
accessible to the experiment. For Atot 27r (50% 
of full sky coverage), one estimates a square win- 
dow of side 4547//V, i.e. 45° for 100 events, 
14° for 1000 events. Since the former number is 
of the order of present world statistics, and the 
latter is the achievement expected by Auger in a 
few years of operations, a binning in windows of 
size 15° represents quite a reasonable choice for 
our forecast. This choice is also suggested by the 
typical size of the observable structures, as can be 
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Figure 1. Galactic coordinate reference frame 
and contours enclosing 68%, 95% and 99% of the 
Auger exposure function, with the corresponding 
declinations. The celestial equator ((5 = 0°) and 
south pole [6 ~ 0, —90°) are also shown. 



seen in the next section. Notice that the GMF, 
that induces at these energies typical deflections 
of about 4° [15], can be safely neglected for this 
kind of analysis. The same remark holds for the 
angular resolution of the experiment. 

Moreover, for a specific experimental set-up 
one must include the proper exposure Wexp, that 
has to be convolved with the previously found 
/lss- We properly parameterize and take into 
account the Auger exposure a;exp(^) as function 
of declination 5 [14,16]: contour plots in galactic 
coordinates are shown in Fig. 1. 

For a given experiment and catalogue, the null 
hypothesis we want to test is that the events 
observed are sampled — apart from a trivial ge- 
ometrical factor — according to the distribution 
/lss '^exp Since we are performing a forecast 
analysis, we will consider test realizations of N 
events sampled according to a random distribu- 
tion on the (accessible) sphere, i.e. according to 
f-^oxp and determine the confidence level (C.L.) 
with which the hypothesis is rejected as a func- 
tion of N . For each realization of N events we 
calculate the two functions 

M 
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{o^ ~_e4j\so\r_ 

[/iso] 

{oi - e4/Lss])' 
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(2) 



(3) 



^the normalization factor is fixed consistently with equa- 
tion (1) 



where Oi is the number of "random" counts in 
the i-th bin f2i, and e,;[/LSs] and e,;[/iso] are the 
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theoretically expected number of events in Oj re- 
spectively for the LSS and isotropic distribution. 
The mock data set is then sampled Af times in 
order to establish empirically the distributions of 
'^LSS "^iso' resulting distribution is 

studied as function of N (plus eventually s, Ecut, 
etc.). 

3. Results 

In Fig. 2 we plot the smoothed maps in 
galactic coordinates of the expected integrated 
flux of UHECRs above the energy threshold 
-Ecut = 3, 5, 7, 9 X 10^^ eV and for slope param- 
eter s = 2.0; the isotropic part has been taken 
into account and the ratio of the isotropic to 
anisotropic part Wiso/ X^fc /Lss(fc) is respectively 
83%,3.6%,<C1%,< 1%. 

Only for Ecut = 3 x 10^^ eV the isotropic 
background constitutes then a relevant fraction, 
since the GZK suppression of far sources is not 
yet present. For the case of interest E'cut = 
5 X 10^^ eV the contribution of Wiso is almost 
negligible, while it practically disappears for 
-Ecut <^ 7 X 10^^ eV. Varying the slope for s = 
1.5,2.0,2.5,3.0 while keeping E^nt = 5 x lO^^eV 
fixed produces respectively the relative weights 
8.0%, 3.6%, 1.8%, 0.9%, so that only for very hard 
spectra Wiso would play a non-negligible role. 

Due to the GZK-effect, as it was expected, the 
nearest structures are also the most prominent 
features in the maps. The most relevant struc- 
ture present in every slide is the Local Superclus- 
ter. It extends along / ~ 140° and / ~ 300° and 
includes the Virgo cluster at / = 284°, b = +75° 
and the Ursa Major cloud at I = 145°, b = +65°, 
both located at z ~ 0.01. The lack of structures 
at latitudes from Z !^ 0° to Z ~ 120° corresponds 
to the Local Void. At higher redshifts the main 
contributions come from the Perseus- Pisces su- 
percluster {I = 160°, 6 = -20°) and the Pavo- 
Indus superclustcr (/ = 340°, 6 = -40°), both at 
z ~ 0.02, and the very massive Shapley Concen- 
tration (Z = 250°, 6 = -1-20°) at 2 ~ 0.05. For a 
more detailed list of features in the map, see the 
key in Fig. 3. 

The EJcut-dependence is clearly evident in the 
maps: as expected, increasing Ecut results in a 




Figure 2. Equal area Hammer-Aitoff projections 
of the smoothed UHECRs arrival directions dis- 
tribution (Eq. (1)) in galactic coordinates ob- 
tained for fixed s = 2.0 and, from the upper to 
the lower panel, for Ecut = 3, 5, 7, 9 x 10^^ eV. 
The smoothing angle is cr = 3°. The contours en- 
close 95%, 68%, 38%, 20% of the corresponding 
distribution. The dark central band corresponds 
to the avoidance zone of the galactic plane. 
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Figure 3. Detailed key of the structures visi- 
ble in the UHECR maps; arbitrary contour lev- 
els. Labels correspond to: (1) Southern exten- 
sion of Virgo and Local Supercluster; (2)Fornax- 
Eridani Cluster; (3) Cassiopea Cluster; (4) Pup- 
pis Cluster; (5) Ursa Major Cloud; (6-7) Pavo- 
Indus and "Great Attractor" region; (8) Centau- 
rus Super- Cluster; (9) Hydra Super-Cluster; (10) 
Perseus Super-Cluster; (11) Abell 569; (12) Pe- 
gasus Cluster; (13-17) Pisces Cluster; (14) Abell 
634; (15) Coma Cluster; (16-18) Hercules Su- 
percluster; (19) Leo Supercluster; (20) Columba 
Cluster; (21) Cetus Cluster; (22) Shapley Con- 
centration; (23) Ursa Major Supercluster; (24) 
Sculptor Supercluster; (25) Bootes Supercluster. 



map that closely reflects the very local universe 
(up to z ^ 0.03 — 0.04) and its large anisotropy; 
conversely, for i?cut — 3, 4 x 10^^ eV, the resulting 
flux is quite isotropic and the structures emerge 
as fluctuations from a background, since the GZK 
suppression is not yet effective. This can be seen 
also comparing the near structures with the most 
distant ones in the catalogue: while the Local 
Supercluster is well visible in all slides, the sig- 
nal from the Perseus-Pisces super-cluster and the 
Shapley concentration is of comparable intensity 
only in the two top panels, while becoming highly 
attenuated for Ecut = 7 x 10"'^^ eV, and almost 
vanishing for Ecut = 9 x 10^^ eV. A similar trend 
is observed for increasing s at fixed i?cut, though 
the dependence is almost one order of magnitude 
weaker. 

Looking at the contour levels in the maps we 
can have a precise idea of the absolute intensity of 
the "fluctuations" induced by the LSS; in partic- 
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Figure 4. The distributions of the estimators 
X^^c, A'Lgg for the cases s ~ 2.0 and for 
N = 200, 1000 events. The distribution arc the 
results of 10.000 monte-carlo simulation like de- 
scribed in the text. 

ular, for the case of interest of Scut = 5 x 10-*^^ eV 
the structures emerge only at the level of 20%- 
30% of the total flux, the 68% of the flux ac- 
tually enclosing almost all the sky. For Ecut = 
7, 9 X 10"'^^ eV, on the contrary, the local struc- 
tures are significantly more pronounced, but in 
this case we have to face with the low statistics 
available at this energies. Then in a low-statistics 
regime it's not an easy task to disentangle the LSS 
and the isotropic distributions. 

In Figure 4 we show the distributions of the 
functions X-^^^ and -^lss introduced in the previ- 
ous section for s = 2.0 and N = 200, 1000. It 
is clear that a few hundreds events are hardly 
enough to reliably distinguish the two models, 
while TV = 800-1000 should be more than enough 
to reject the hypothesis at 2-3 cr, independently of 
the injection spectrum. Steeper spectra however 
slightly reduce the number of events needed for a 
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Figure 5. The first multipoles of the dimension- 
less angular power spectrum of the 5 x 10^^ eV 
map from the fitting procedure; also showed, by 
comparison, the cosmic variance errors slightly 
shifted for clarity. 

given C.L. discrimination. Respect to the avail- 
able data our results clearly show that AGASA 
statistics (only 32 data at > 5 x 10^^ eV in 
the published data set [17], some of which falling 
inside the mask) is too limited to draw any firm 
conclusion on the hypothesis considered. 

4. Angular Power Spectrum 

The previous forecast analysis is completely 
independent from an harmonic analysis of the 
anisotropics; however, a multipole decomposi- 
tion and the related power spectrum calculation 
greatly helps in the understanding of the UHE- 
CRs anisotropics and, indeed, a dipole anisotropy 
study has been the aim of various works [18-20]; 
we cover the gap in this section. 

In the following we will concentrate on the 
-E-cut = 5 X 10^^ eV map; moreover. Auger is not 
expected to be very sensitive to multipoles higher 
than 10, even after ten years of operations, so we 
limit the analysis to the first few I. 

The main problem in calculating the Spherical 
Harmonics expansion of the UHECRs map is the 
presence of the galactic avoidance zone; however, 
in the first few Vs the typical angular scale as- 
sociated to the given multipole (of the order of 
~ tt/Z) is larger of the average angular extension 
of the galactic cut so that an extrapolation in this 



zone is justified. We implement this method fit- 
ting to the UHECRs map the sum of the first 
harmonics till Zmax- The corresponding recon- 
structed map adding up to the first ^=10 and 1=2 
harmonics is showed in Fig. 6 together with the 
original masked map. 

Of course, filling the mask requires an extrapo- 
lation from the near regions and this is somewhat 
questionable in relation to the fact that some un- 
known cluster can be hidden behind our galaxy 
changing in principle all the multipoles pattern; 
a not well quantifiable error is indeed related to 
this ignorance though some more sophisticated 
astronomical analysis indicate that it is unlikely 
that relevant structures hide behind the galactic 
plane [21-23] so that at this angular resolution 
{I = 10 corresponds to 180/10 w 20° resolution) 
an extrapolation is quite justified. 

In this respect the situation is quite different 
from the CMB case where the anisotropics are 
isotropic and gaussian so that even a limited re- 
gion of sky represents a fair statistical sample of 
the whole CMB sky, while the various multipoles 
are independent realization of the underlying ran- 
dom field. The UHECRs field is instead related 
to local and non-linear structures and it is then 
far from being gaussian; moreover the presence of 
possible unknown clusters in the avoidance zone 
introduces a correlation among the various mu- 
tipoles so that the naive error bars derived from 
the fit are not really trustable. We instead quan- 
tify the error looking at the spread of the best fit 
values obtained varying the degrees of freedom 
involved into the fit, related in its turn to the 
chosen ^max- We varied Imax in the range 1-12 
that we found to be optimal. 

The resulting angular power spectrum for the 
first 10 multipoles of the adimensional map 
fhssi^, Ecut) / </lss> is shown in Fig. 5; note 
that the true errors are not affected by the cosmic 
variance given that we know the complete pattern 
of the angular anisotropics and not only their sta- 
tistical properties. The hypothetical errors from 
the cosmic variance are, however, showed by com- 
parison: as expected the true errors on the very 
wide structures of the dipole and quadrupole are 
quite small compared to cosmic variance; the er- 
rors, then, rapidly degrade when the scale related 
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Figure 6. (a) The starting UHECRs map and 
mask in equatorial coordinates, (b) The recon- 
structed map from the first 10 multipoles. (c) As 
above for dipole-|-quadrupole. 



to the given I approaches the size of the galactic 
cut so that already for Z = 8 — 10 the errors are 
almost comparable to the cosmic variance. As 
expected the dipolc gives a prominent contribu- 
tion though a similar power is predicted for the 
quadrupole; this is easily explained in terms of the 
original map where besides the Super-Galactic 
plane there is an important contribution from 
the P-P cluster, which induces a mixed dipolar- 
quadrupolar pattern. This contribution is almost 
an order of magnitude greater respect to the other 
multipoles so that, likely, it is also the dominant 
contribution in assessing the minimum statistics 
needed to detect the anisotropics like determined 
in the previous paragraphs. Since the dipole and 
quadrupole are the easiest anisotropies to look for 



Table 1 

Low I multipole moments in equatorial coordi- 
nates. a;_m moments are obtained through the 
reality condition ai-m = Cij[m- 

in UHECR analyses, this result suggests that the 
1=1+1=2 pattern should slowly emerge from the 
cosmic rays data as the first detectable large scale 
anisotropy. Given the relevant role represented 
by the quadrupole and dipole we give in table 1 
the expression of the aim and the related errors, 
in equatorial coordinates. Finally, it is worth to 
notice the scale of the multipole anisotropies of 
the order ~ 0.06 that translates into an overall 
anisotropy of the order \/0.06 « 25% in good 
agreement with what found with analysis of the 
full map. 

As check of the method we perform an inde- 
pendent estimate of the power spectrum using 
the Master algorithm [24] , a method typically em- 
ployed in similar CMB analysis. This method is 
well suited for small scale anisotropies and is quite 
complementary to the fitting procedure; we ver- 
ified that the two spectra superimpose and join 
smoothly making us confident on the reliability 
of the result. 

5. Summary and conclusion 

UHECRs offer the fascinating possibility to 
make astronomy with high energy charged parti- 
cles; in particular, besides small scale clustering, a 
large scale anisotropy is expected correlating with 
the local cosmological LSS. However, we showed 
with a approach that a huge statistics of sev- 
eral hundreds data is required to start testing the 
model at Auger South. 

Untill now, the lack of UHECR statistics has 
seriously limited the usefulness of such a kind 
of analysis. However, progresses are expected in 
forthcoming years. The Southern site of Auger 
is already taking data and, once completed, the 
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total area covered will be of 3000 km^, thus im- 
proving by one order of magnitude present statis- 
tics in a couple of years [25] . The idea to build a 
Northern Auger site strongly depends on the pos- 
sibility to perform UHECR astronomy, for which 
full sky coverage is of primary importance. To 
this aim, the Japanese- American Telescope Array 
in the desert of Utah is expected to become op- 
erational by 2007 [26] offering almost an order of 
magnitude larger aperture per year than AGASA 
in the Northern sky. 
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